A mathematical model for optimizing the indications of liver transplantation in patients with hepatocellular carcinoma

Abstract Background The criteria for organ sharing has developed a system that prioritizes liver transplantation (LT) for patients with hepatocellular carcinoma (HCC) who have the highest risk of wait-list mortality. In some countries this model allows patients only within the Milan Criteria (MC, defined by the presence of a single nodule up to 5 cm, up to three nodules none larger than 3 cm, with no evidence of extrahepatic spread or macrovascular invasion) to be evaluated for liver transplantation. This police implies that some patients with HCC slightly more advanced than those allowed by the current strict selection criteria will be excluded, even though LT for these patients might be associated with acceptable long-term outcomes. Methods We propose a mathematical approach to study the consequences of relaxing the MC for patients with HCC that do not comply with the current rules for inclusion in the transplantation candidate list. We consider overall 5-years survival rates compatible with the ones reported in the literature. We calculate the best strategy that would minimize the total mortality of the affected population, that is, the total number of people in both groups of HCC patients that die after 5 years of the implementation of the strategy, either by post-transplantation death or by death due to the basic HCC. We illustrate the above analysis with a simulation of a theoretical population of 1,500 HCC patients with tumor size exponentially. The parameter λ obtained from the literature was equal to 0.3. As the total number of patients in these real samples was 327 patients, this implied in an average size of 3.3 cm and a 95% confidence interval of [2.9; 3.7]. The total number of available livers to be grafted was assumed to be 500. Results With 1500 patients in the waiting list and 500 grafts available we simulated the total number of deaths in both transplanted and non-transplanted HCC patients after 5 years as a function of the tumor size of transplanted patients. The total number of deaths drops down monotonically with tumor size, reaching a minimum at size equals to 7 cm, increasing from thereafter. With tumor size equals to 10 cm the total mortality is equal to the 5 cm threshold of the Milan criteria. Conclusion We concluded that it is possible to include patients with tumor size up to 10 cm without increasing the total mortality of this population.

Mechanical Alterations in airway smooth muscle after interaction with indoor pollutant – A mathematical model.

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully understood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholine (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle ¯bers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, incorporates the phenomenology of both MCh and FA and reproduces experimental results observed with in vitro exposure of smooth muscle to FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells in a tissue level model. The model can also be used in different biological scales.

Mechanical alterations in airway smooth muscle after interaction with indoor pollutant - A mathematical model

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully w1derstood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholinc (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle fibers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, ir1corporates the phenomenology of both MCh and FA and reproduces experirnental results observed with ir1 vitro exposure of smooth muscle to .FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells ir1 a tissue level model. The model can also be used in different biological scales.

Effects of arterial oxygen tension and cardiac output on venous saturation: a mathematical modeling approach

OBJECTIVES: Hemodynamic support is aimed at providing adequate O-2 delivery to the tissues; most interventions target O-2 delivery increase. Mixed venous O-2 saturation is a frequently used parameter to evaluate the adequacy of O-2 delivery. METHODS: We describe a mathematical model to compare the effects of increasing O-2 delivery on venous oxygen saturation through increases in the inspired O-2 fraction versus increases in cardiac output. The model was created based on the lungs, which were divided into shunted and non-shunted areas, and on seven peripheral compartments, each with normal values of perfusion, optimal oxygen consumption, and critical O-2 extraction rate. O-2 delivery was increased by changing the inspired fraction of oxygen from 0.21 to 1.0 in steps of 0.1 under conditions of low (2.0 L.min(-1)) or normal (6.5 L.min(-1)) cardiac output. The same O-2 delivery values were also obtained by maintaining a fixed O-2 inspired fraction value of 0.21 while changing cardiac output. RESULTS: Venous oxygen saturation was higher when produced through increases in inspired O-2 fraction versus increases in cardiac output, even at the same O-2 delivery and consumption values. Specifically, at high inspired O-2 fractions, the measured O-2 saturation values failed to detect conditions of low oxygen supply. CONCLUSIONS: The mode of O-2 delivery optimization, specifically increases in the fraction of inspired oxygen versus increases in cardiac output, can compromise the capability of the "venous O-2 saturation" parameter to measure the adequacy of oxygen supply. Consequently, venous saturation at high inspired O-2 fractions should be interpreted with caution.

A heat transfer model of the human upper limbs

A steady state multi-segmented heat transfer model of the human upper limbs was developed. The main purpose was to evaluate the impact of blood flow through superficial veins and subcutaneous vascular structures in the palm of the hands over the heat transfer between the limbs and the environment. The distinguishing feature of the model is the inclusion of a detailed circulatory network composed of vessels with diameter larger than 1 mm. The model was validated by comparing its results from exposures to a hot, a neutral, and a cold environment to experimental data presented in the literature. (C) 2011 Elsevier Ltd. All rights reserved.

Forced Desynchronization of Activity Rhythms in a Model of Chronic Jet Lag in Mice

We studied locomotor activity rhythms of C57/Bl6 mice under a chronic jet lag (CJL) protocol (ChrA(6/2)), which consisted of 6-hour phase advances of the light-dark schedule (LD) every 2 days. Through periodogram analysis, we found 2 components of the activity rhythm: a short-period component (21.01 +/- 0.04 h) that was entrained by the LD schedule and a long-period component (24.68 +/- 0.26 h). We developed a mathematical model comprising 2 coupled circadian oscillators that was tested experimentally with different CJL schedules. Our simulations suggested that under CJL, the system behaves as if it were under a zeitgeber with a period determined by (24 -[phase shift size/days between shifts]). Desynchronization within the system arises according to whether this effective zeitgeber is inside or outside the range of entrainment of the oscillators. In this sense, ChrA(6/2) is interpreted as a (24 - 6/2 = 21 h) zeitgeber, and simulations predicted the behavior of mice under other CJL schedules with an effective 21-hour zeitgeber. Animals studied under an asymmetric T = 21 h zeitgeber (carried out by a 3-hour shortening of every dark phase) showed 2 activity components as observed under ChrA(6/2): an entrained short-period (21.01 +/- 0.03 h) and a long-period component (23.93 +/- 0.31 h). Internal desynchronization was lost when mice were subjected to 9-hour advances every 3 days, a possibility also contemplated by the simulations. Simulations also predicted that desynchronization should be less prevalent under delaying than under advancing CJL. Indeed, most mice subjected to 6-hour delay shifts every 2 days (an effective 27-hour zeitgeber) displayed a single entrained activity component (26.92 +/- 0.11 h). Our results demonstrate that the disruption provoked by CJL schedules is not dependent on the phase-shift magnitude or the frequency of the shifts separately but on the combination of both, through its ratio and additionally on their absolute values. In this study, we present a novel model of forced desynchronization in mice under a specific CJL schedule; in addition, our model provides theoretical tools for the evaluation of circadian disruption under CJL conditions that are currently used in circadian research.

An in-depth analysis of the HIV-1/AIDS dynamics by comprehensive mathematical modeling

This work presents major results from a novel dynamic model intended to deterministically represent the complex relation between HIV-1 and the human immune system. The novel structure of the model extends previous work by representing different host anatomic compartments under a more in-depth cellular and molecular immunological phenomenology. Recently identified mechanisms related to HIV-1 infection as well as other well known relevant mechanisms typically ignored in mathematical models of HIV-1 pathogenesis and immunology, such as cell-cell transmission, are also addressed. (C) 2011 Elsevier Ltd. All rights reserved.

Numerical simulation of an ethanol turbulent spray flame with RANS and diffusion combustion model

A detailed numerical simulation of ethanol turbulent spray combustion on a rounded jet flame is pre- sented in this article. The focus is to propose a robust mathematical model with relatively low complexity sub- models to reproduce the main characteristics of the cou- pling between both phases, such as the turbulence modulation, turbulent droplets dissipation, and evaporative cooling effect. A RANS turbulent model is implemented. Special features of the model include an Eulerian– Lagrangian procedure under a fully two-way coupling and a modified flame sheet model with a joint mixture fraction– enthalpy b -PDF. Reasonable agreement between measured and computed mean profiles of temperature of the gas phase and droplet size distributions is achieved. Deviations found between measured and predicted mean velocity profiles are attributed to the turbulent combustion modeling adopted